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ANALYSIS AND CORRELATION OF HEAT-TRANSFER 
COEFFICIENT AND FRICTION FACTOR DATA 

FOR DILUTE GAS-SOLID SUSPENSIONS 

by Robert Pfeffer,* Salvatore Rossetti, t and Seymour Lieblein 

Lewis Research Center 

SUMMARY 

The available literature on both the heat transfer and the pressure drop associated 
with the flow of a dilute gas -solid suspension was reviewed and analyzed. This study 
obtained workable correlations for predicting some of the basic thermophysical proper- 
ties and both the heat-transfer coefficient and the friction factor for the flow of suspen- 
sions in tubes. The basic properties of suspensions investigated included bulk density, 
specific heat, thermal conductivity, viscosity, specific heat ratio, and sonic velocity. 

The variables considered included the loading ratio (pounds of solid per pound of gas), 
the ratio of particle to gas specific heat, the gas Reynolds number, and the particle 
diameter. 

Two possible correlations for the convective heat-transfer coefficient of a suspen- 
sion resulted from this analysis. Both correlations indicated that the ratio of suspension 
to pure gas heat-transfer coefficient increased with loading ratio and specific heat ratio, 
decreased with increasing gas Reynolds number, and was essentially unaffected by 
particle diameter for the range of conditions considered. 

*Professor of Chemical Engineering, City College of the City University of New 
York. 

^Graduate Student, City College of the City University of New York. 
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The Reynolds analogy was applied to these heat -transfer relations, and a correlation 
for the ratio of suspension to gas equivalent friction factor was obtained. The ratio of 
the suspension friction factor to the pure gas friction factor was found to increase with 
loading ratio and decrease with gas Reynolds number and to be essentially unaffected by 
particle diameter. This correlation was found to be in good agreement with both the ex- 
perimental data and an analytical study based on an eddy viscosity model for a suspension. 


INTRODUCTION 

Space power generation systems must be capable of continuously generating power 
for long periods of time. One power generating system currently under consideration is 
the indirect-conversion closed-loop heat engine. In this system heat is generated in a 
nuclear or solar source and rejected by a radiator, with power obtained from a turbine. 
However, because of mission limitations the powerplant specific weight must be kept low. 

The Brayton cycle using a pure inert gas as the working fluid has been considered as 
a potential source of space power (refs. 1, 2, and 3). The basic technology available 
from developments in the gas turbine and power production fields can be advantageously 
utilized in Brayton cycle space power applications. However, because of the low temper- 
ature levels in the radiator the Brayton cycle waste heat rejection unit becomes very large 
especially when a gas radiator is used in the power loop. It is apparent, therefore, that 
if the area and weight of the radiator could be reduced, this system would become more 
promising for future space missions. 

In one approach to circumvent the effects of the inherently low heat -transfer coeffi- 
cients in radiator gas flow, a second loop containing a liquid-flow radiator in conjunction 
with a gas-to-liquid heat exchanger and circulating pump can be used instead of the gas 
flow radiator. However, such an arrangement involves increased complexity and an addi- 
tional power drain. Another method of decreasing the radiator area, and hence its weight, 
is to increase the heat-transfer coefficient of the gas working fluid to allow a greater heat 
transfer per unit of radiation area. Recent studies have shown that the addition of solid 
particles to a turbulently flowing gas will increase the rate of heat transfer between the 
gas and its surroundings (refs. 4 to 6). This increase in heat-transfer rate results from 
two separate effects: an increase in the volumetric heat capacity of the working fluid (be- 
cause the axial temperature difference is decreased), and an increase in the gas side heat- 
transfer coefficient. Because of these desirable properties, the use of a gas-solid sus- 
pension of small particles is being considered as a coolant for nuclear reactors, and 
might also be useful as a working fluid in a Brayton cycle space power generating system. 

In these applications of dilute flowing gas-solid suspensions, it is necessary to be 
able to predict the suspension heat-transfer coefficient and pressure drop so that the heat 
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rejection and pumping equipment can be adequately designed. However, readily usable 
correlations for predicting both heat-transfer rates and pressure drop of flowing suspen- 
sions are not available. In an attempt to establish such correlations, a review and ana- 
lysis of available experimental data, correlations, and theory pertaining to heat -transfer 
coefficients and pressure drops of gas -solid suspensions was undertaken. Some of the 
related thermodynamic and transport properties of gas-solid suspensions were also re- 
viewed. The analysis was limited to dilute suspensions of small particles (loading ratios 
less than 10 and particle diameters less than 150/i.) since these are of primary interest 
in application to the Brayton cycle. 


THERMOPHYSICAL PROPERTIES 

The performance of a flowing dilute gas -solid suspension (heat transfer, pressure 
drop, etc. ) will be dependent on the fundamental thermophysical properties of the suspen- 
sion. In general, these properties of the suspension will be different from those of the 
pure carrier gas. Thus, before proceeding with the analysis of heat-transfer and pres- 
sure drop characteristics, it is necessary to define the principal thermophysical proper- 
ties of the suspension. Some of the more fundamental thermodynamic and transport prop- 
erties of interest in a gas-solid suspension are density, specific heat, thermal conduc- 
tivity, viscosity, specific heat ratio, and sonic velocity. In many cases, the properties 
of the suspension will be related to those of the pure gas and will invariably be a function 
of the loading ratio, defined as the ratio of the mass flow rate of solid to the mass flow 
rate of the pure gas. 


Bulk Density 

The bulk density of the suspension p„ and the densities of the pure substances 
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(particles and gas) are related to the loading ratio r\ by the equation 


where 



77 + 1 

i?(pg/P p ) + 1 



(la) 


(lb) 


I' 


3 



Since Pg/P p is generally small /Pg/P p < 10 and if rj <10, then f?(Pg/p p ) « 1, and 
equation (la) reduces to \ ' 


1 + TJ 


(lc) 


(Symbols are defined in the appendix. ) 


Specific Heat 


If a weighted average for the specific heat of the suspension is assumed, 

( c p) g 1 + * 

Since by definition 



Then 


^ A _ 1 + 5r? 


( 2 ) 


(3) 


A plot of against r] for a suspension of graphite particles in helium, neon, 

and argon is given in figure 1. At large values of r], the heat capacity ratio 
asymptotically approaches the value of 6 for each suspension. 


( c p) s /H 


Thermal Conductivity 

Gorring and Churchill (ref. 7) have presented an analytical equation that agrees quite 
well with experimentally measured thermal conductivities of dispersions. They have 
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shown that for ij < 20 the effective thermal conductivity of a suspension is given by 

k s _ j3 - € 

k B + e 
g p 


O Q 

where e is the fractional solid volume (ft solid/ft suspension) given by 

vp c 


e = 


g 


n 




and 



Hence 





Figure 2. - Effect of loading ratio on ratio of thermal conductivity of suspension to thermal 
conductivity of pure gas for systems of graphite suspended in argon, neon, and helium. 
Temperature, 700° R; density of graphite, 105 pounds per cubic foot; thermal conductivity 
of graphite, 100 Btu/(hr)(ftM°R). 



( 6 ) 


Results of sample calculations for k g /kg for several gases with graphite particles 
are shown in figure 2. For t] < 10, the conductivity of the suspension was not seen to be 
very different from that of the pure gas for the cases considered. 


Viscosity 

Very little experimental or analytical work has been done in determining the viscosity 
of a gas-solid suspension. The only reference found for experimental work on the viscos- 
ity of gas-solid suspensions was done by Sproull (ref. 8) at low loading ratios fa < 1) with 
dusty air. The results that Sproull found at these low loading ratios indicate that the 
viscosity of a gas -solid suspension is actually lower than that of a pure gas; however, at 
these very low particle concentrations the concept of the viscosity of a suspension be- 
comes nebulous. 

In an analytical study Happel (ref. 9), by omitting the inertia terms in the steady - 
state Navier -Stokes equations, derived the following relation between relative viscosity 
and solid concentration that is in good agreement with existing data at high loading ratios 

( V > 20 ): 


6 




Figure 3. - Effect of loading ratio on ratio of viscosity of suspension to viscosity of pure gas 
for systems of graphite suspended in argon, neon, and helium. Temperature, 700° R; 
density of graphite, 105 pounds per cubic foot 


fl 

— =1 + 5.5 (7a) 

\i r _ 

S 

where e is the fraction of volume occupied by solids as given by equation (5), and \f/ is 
a factor characteristic of the loading ratio. 

For 77 < 10, x(/ — 1. 0, and equation (7a) reduces to 

M 

= 1 + 5. 5 e (7b) 

which is identical in form to the Einstein equation (ref. 10), except for the difference in 

the constant, 5. 5 for equation (7b) as compared with 2. 5 for the Einstein equation. 

Sample results for jtz_/ ju_, as given by equation (7b), are shown in figure 3 for three 

s g 

gases with graphite particles. As can be seen in the figure, the suspension viscosity may 
be substantially different from the gas viscosity for heavy gases at high pressure. 

Specific Heat Ratio 

The effective isentropic specific heat ratio v for a suspension can be readily obtained 

s 

from an energy balance if the temperature and velocity of the particles in the suspension are 
assumed to be the same as the temperature and velocity of the gas (refs. Hand 12). Based 
on these assumptions and using the ratio of the specific heat of the particles to the specific 
heat of the gas 6 , the ratio of specific heat ratios in terms of the loading ratio 77 becomes 
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sions. 

Is. = 1 + ft? (8) 

y g 1 + y g 571 

For large loading ratios, 7] > 10, y g — 1. 0, and the suspension behavior approaches that 
of an incompressible fluid. Figure 4 is a plot of the ratio of the isentropic specific heat 
ratios y g /y as a function of v for suspensions of graphite particles in inert gases. 

The figure indicates that for neon and argon y g /y g decreases rapidly as t? increases. 
For the graphite -helium system, the decline in y g /y g with r\ is less pronounced be- 
cause of the low value of 6 for this system. For all three systems y g /y g approaches 
the asymptotic value of 0. 6. 

Since solid particles can do no expansion work, their presence in components in- 
volving flow expansion, such as rocket nozzles or gas turbines, can only decrease the 
effectiveness of the expansion process in converting thermal energy to kinetic energy. 

The particles are accelerated entirely by drag forces associated with lag or slippage of 
the particles relative to the expanding gas; therefore, some performance loss relative to 
the calculated ideal no-slip expansion must be present when dealing with solid particles. 
Experience has shown that both velocity lag and thermal lag become significant as the 
size of the particles and the particle loading ratio increase; therefore, the value of y_ 
given by equation (8) is at best an approximation or a limiting value. 

A simplified one -dimensional analysis by Kleigel (ref. 13) for a gas-solid suspension 
flow in a rocket nozzle shows that y_ increased from the equilibrium value given by 

O 

equation (8) as the particle was accelerated but was rather insensitive to particle lag for 
small lags. This analysis indicates that the use of equation (8) to estimate y g may be a 
reasonable approximation even for nonequilibrium flow when the particle size is small. 

In a nozzle expansion, for example, 1- to 2-micron particles were found to follow the gas 
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closely, whereas 10-micron particles had a significant lag (ref. 13) in the throat of the 
nozzle. 


Sonic Velocity 


Assuming that Vp = Vg 
as (ref. 12) 


and T = T 
P g 


gives the velocity of sound in the suspension 


“s/vnA 178 

a g \ 1+ ” / 


( 9 ) 


where a is the velocity of sound in the pure gas at the same temperature. Equation (9) 
shows that the sonic speed of the suspension a? g is always less than that of a pure gas. 

A plot of a g / Qfg against r\ is given in figure 5. The figure shows that the sonic veloc- 
ity of the suspension, as given by equation (9), decreases rapidly as 7j is increased; 
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however, how this variation will affect the compressibility and choking of the flow of a 
dilute suspension in which only the gas carrier is essentially the compressible phase is 
not clear. 


CONVECTIVE HEAT TRANSFER 

There appear to be two possible mechanisms explaining why the addition of solid 
particles to a turbulently flowing gas will increase the rate of heat transfer between the 
gas and its surroundings. In the first mechanism this increase in heat-transfer rate is 
attributed to the particles eroding away some of the stagnant gas film at the containing 
surface and thereby reducing the resistance to heat transfer between the particles and the 
surface. In the other suggested mechanism the heat transfer from the particles to a sur- 
face is mainly affected by the mechanical contact with the wall of particles or groups of 
particles that are continuously replaced at the exchanger surface. In this model, parti- 
cles at the bulk temperature of the suspension are swept to the heat-transfer surface by 
the turbulent motion of the suspension. Here, the relatively high temperature gradient 
between the heat-transfer surface and the particle adjacent to it transfers heat by con- 
duction at the surface -particle contact. As the temperature of the particle begins to 
approach that of the heat-transfer surface, the driving potential is reduced and heat flows 
less rapidly. Thus, achievement of a large heat-transfer coefficient is dependent on 
vigorous agitation within the bulk of the suspension and on short residence times for the 
particle at the heat-transfer surface. 

If the increased heat transfer is due only to the mechanical shear of the laminar gas 
film by the particles, the coefficient might appear to depend on the loading ratio but 
should be independent of the specific heat and thermal conductivity of the solid. On the 
other hand, if the increase in heat-transfer rate is due to contact of particles with the 
wall surface, then it would be expected that the thermal conductivities of the solid and gas 
play a significant role in determining the heat-transfer coefficient. Further experimental 
studies on a number of dissimilar solids and gases might elucidate the true mechanism of 
the action of the solids. 

Since there is no general agreement on the mechanism by which heat-transfer coef- 
ficients of gases are increased by the addition of solids, no attempt will be made to 
validate either concept. Rather, it will be the purpose of this section to review the ana- 
lytical and experimental literature dealing with heat-transfer coefficients of gas -solid 
suspensions in order to obtain a working correlation for the ratio of the heat-transfer 
coefficient of the suspension to the heat-transfer coefficient of the pure gas h g /hg for 
fully developed turbulent flow in circular tubes. 

The ratio of suspension to pure gas heat-transfer coefficient can be obtained by 
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relating the measured coefficient at various loading ratios to the measured coefficient at 
rj = 0 or by relating the suspension Nusselt number to the pure gas Nusselt number. 

The Nusselt number of the suspension is defined as 


Nu 

s 



( 10 ) 


For the pure gas (ref. 14) 


Nu = = 0. 023 Re 8 ' 8 Pr 8 ' 3 

g k g g 

K g 


for cooling 


and 


h D 

NUg = -S- = 0. 023 Re 8, 8 Pr 8 " ^ for heating 
k g 


In terms of Nusselt numbers 


(ID 


( 12 ) 



Nu k 
s s 

Nu„ k 
g g 


(13a) 


However, for the test conditions of the available experimental data and correlations, 

k /k , as computed by equation (6), was generally sufficiently close to unity (<1. 05) 
s g 

for rj < 10 so that for simplicity 


h 

g 


Nu 

s 

N u 

g 


(13b) 


Note in equation (13b) that the Nusselt number ratio should be evaluated for the same gas 
conditions for which the heat-transfer coefficient ratio is desired (e.g. , same gas pres- 
sure, temperature, and volumetric flow rate for the suspension and pure gas flows). 


Data Sources 

Many experimental studies have attempted to investigate the effect of different 
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TABLE I. - SUMMARY OF OPERATING CONDITIONS UTILIZED IN HEAT-TRANSFER WORK CITED 


Particle 

Gas 

Specific 

heat 

ratio, 

6 

Particle 

diameter, 

d p> 

V 

Loading 

ratio, 

V 

Gas 

velocity, 

V 

ft/sec 

Tube inside 
diameter, 
D, 
in. 

Gas Reynolds 
number, 

Re g 

Ratio of 
tube length 
to tube 
diameter, 
L/D 

Test 

condition 

Fluid 

temperature, 

T, 

°F 

Source 

Silica- 

alumina 

catalyst 

Air 

a l. 15 

10 to 210 
(50{i avg) 

0 to 13. 3 

38. 5 to 81 

0.7 

13 000 to 
28 000 

46.6 

Heating 

75 to 102 

Farbar and 
Morley 
(ref. 4) 

Glass 

Air 

0.797 

30, 70, 
140, 200 

0 to 10. 5 

41 to 71 

0.67 

15 300 to 
26 500 

47 

Heating 


Farbar and 
Depew 
(ref. 15) 


Silica - 

alumina 

catalyst 

Air 

0.84 

10 to 210 
(50 p. avg) 

2 to 446 

1.4 to 
61.5 

0. 688 to 
1.5 

178 to 
25 400 

240 

Cooling 

1050 

Danziger 
(ref. 5) 

Graphite 

Nitrogen, 

helium, 

carbon 

tetra- 

fluoride, 

carbon 

dioxide 

0.315 to 
1.50 

1 to 5 

0 to 90 

20 to 
200 

0.313 to 
0.875 

2000 to 
400 000 

100 to 
289 

Heating 

90 to 1100 

Schluderberg, 
Whitelaw and 
Carlson (ref. 18), 
and Franklin Inst, 
(ref. 19) 

Graphite 

Air 

0. 84 to 
1.03 

150, 1440, 
400, 770 
1160,2080 

0 to 32 

25 to 95 

0.473 to 
1.3 

7000 to 
65 000 

50 

Heating 


Gorbis and 
Bakhtiozin 
(ref. 6) 


Glass 

and 

lead 

Air 

0:-80 to 
0. 133 

30, 200 

0 to 3 

40 to 85 

0.71 

15 000 to 
30 000 

46.4 

Heating 

100 

Tien and 
Quan (ref. 23) 

Glass 

Air 

0.80 

40, 284, 
100, 450, 
150, 270 

45 to 
1600 

0. 8 to 15 

1 to 4 

225 to 
16 800 

Fluidized 

bed 

Heating 

120 to 
230 

Mickley and 
Trilling 
(ref. 24) 

Glass 

and 

copper 

Air 

0.80, 

0.675 

50, 534, ' 

224, 586 
392, 356 

0. 1 to 
16.5 

10 to 34 

li to 2 
2 

7250 to 
33 000 

60 to 
80 

Heating 

85 to 
350 

List (ref. 25) 


a This value has been disputed by Danziger (ref. 5). 


physical parameters on the heat-transfer coefficient of gas -solid suspensions. A brief 
summary of the more important and up-to-date research that has been done in this field 
will now be presented. For convenience, a summary of the principal operating conditions 
and physical properties used by each investigator appears in table I. 

Farbar and Morley . - An exploratory paper by Farbar and Morley (ref. 4) studied 
the effect on heat-transfer rate of adding silica -alumina catalyst particles to air flowing 
in a tube 0. 7 inch in diameter. The particle size was not uniform but varied from 
10 to 210 microns with half of the particles smaller than 50 microns. The range of gas 
Reynolds number Re g varied between 13 000 and 28 000, and the loading ratio tj was 
varied between 0 and 13. 3. Farbar and Morley found that for a solid loading ratio of 
unity or less there is practically no effect on the value of the heat-transfer coefficient, 


12 



whereas at i\ greater than 1, the heat-transfer coefficient increased as rj was 
increased. They also found that increasing the gas Reynolds number had a smaller ef- 
fect on the heat-transfer coefficient of the suspension than on that for a pure gas. This 
is revealed in their final correlation, which is given by 

Nu = 0. 14 Re°‘ 6 tj°- 45 (14) 

s g ' 

If equation (14) is divided by equation (12) and the Prandtl number for air is taken as 
0. 75, the heat-transfer coefficient ratio becomes 



6.8 Re" 0,2 fj 0 - 45 


(15) 


The ratio of specific heats 6 and the particle diameter Dp do not appear in their cor- 
relation since all of Farbar and Mor ley’s experiments were run with only one gas and 
one suspension material. 

Farbar and Depew. - In a paper by Farbar and Depew (ref. 15), the work of Farbar 
and Morley was extended in order to attempt to isolate and determine the effect of parti- 
cle diameter on the heat transfer of gas -solid suspensions. In this study spherical glass 
particles of four uniform sizes (30, 70, 140, and 200fi) were added to air flowing verti- 
cally in a borosilicate glass tube. The loading ratio was varied from 0 to 10, while the 
gas Reynolds number was held constant at 15 300, 19 700, and 26 500. The results in- 
dicated a substantial increase in the gas -side heat-transfer coefficient for the 30-micron 
particles, a moderate increase for the 70-micron particles, a slight increase for the 
140-micron particles, and essentially no increase for the 200-micron particles for the 
range of variables covered. The results of Farbar and Depew were not correlated al- 
though they did clearly indicate a decrease in heat-transfer coefficient with increase in 
particle diameter. The data also indicated a decrease in h g /hg with increasing gas 
Reynolds number at a given loading ratio and particle size. 

Danziger . - W. J. Danziger (ref. 5) presented a report in which a correlation of 
the heat-transfer coefficients of silica -alumina cracking catalyst of about 50-micron 
average particle size in vertical transport was obtained. His data were obtained on re- 
cycle catalyst coolers of two designs, both used vertical, single-tube pass, removable- 
bundle, fire -tube boilers with air-catalyst mixture flowing upward through the tubes. 

The correlation that Danziger obtained based on two sets of data, his own data for cooling 
and the data of Farbar and Morley for heating is 
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( 16 ) 


^= 0.0784Re“' 6 V- 45 

k g 

g 


Converting this correlation to h _/h by dividing by equation (11) and taking Pr = 0. 75 

S g g 

gives 


^s.TRe - 0 '^ 0 - 45 


g 


(17) 


This correlation covers gas Reynolds numbers from 178 to 25 400, loading ratios of 
2 to 446, and tube inside diameters from 0. 689 to 1. 497 inches. The correlation is based 
on both commercial data as well as the laboratory data of Farbar and Morley; however, 
it is limited to one average particle diameter (50/u) and one value of 6. 

Schluderberg, Whitelaw, and Carlson. - In an effort to study the properties of gas- 
eous suspensions as reactor coolants, an extensive research program was conducted that 
is fully described by a series of reports (e. g. , refs. 16 and 17) and is summarized in a 
paper by Schluderberg, Whitelaw, and Carlson (ref. 18). These studies of heat-transfer 
and pressure drop of suspensions of 1- to 5-micron graphite particles in carbon dioxide, 
helium, nitrogen, and carbon tetrafluoride were carried out at gas pressures between 
30 and 130 pounds per square inch gage, gas temperatures between 90° and 1100° F, and 
suspension densities of up to 8 pounds per cubic foot (rj up to 90). 

The heat-transfer data of references 16 to 18 for tube inside diameters from 0. 313 

3 5 

to 0. 875 inch and gas Reynolds numbers from 2x10 to 4x10 were correlated by the 
following equation: 

Nu = 0. 02 Re®' 8 Pr®' 8 (1 + 6*7)®' 45 (18) 

° D 0 

Therefore, by dividing equation (18) by equation (12), the ratio of the heat-transfer 

coefficient of the suspension to the heat-transfer coefficient of the pure gas (assuming 

an average Pr of 0. 75) can be approximated by 
& 


— = 0.78 (1 + 677)®' 45 
h„ 


(19) 


Equation (19) indicates that the increase in heat-transfer coefficient ratio depends pri- 
marily on the specific heat parameter 6 and the loading ratio ij. It shows no effect of 
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gas Reynolds number or particle size on the heat-transfer coefficient ratio. 

Franklin Institute . - The Franklin Institute (ref. 19) conducted a critical evaluation 
of the work reported earlier by Schluderberg, et al. by starting with the original raw 
data and recalculating all of the heat-transfer coefficients and other correlation param- 
eters. They found that the data, although believed generally unreliable because of faulty 
heat balances, could be correlated by any of the following equations: 


NUg 


= 0.0205 (1 + 7f)“ 8 ' 4 



Pr 


0.4 

s 


( 20 ) 


Nu = 0.4 Re 8 * 5 Pr 8 ' 5 (l + fy) 8 ' 45 (21) 

o o 

Nu g = (Re s Pr g )°- 42 (22) 

Nu = 0. 023 (1 + v)~°‘ 7 Re?' 8 Pr?' 4 (23) 

o o S 

and 

Nu s = 0. 015 (1 + t?)-°' 4 Re 8 ' 8 Pr 8 * 4 (24) 


where T g is the bulk temperature of the suspension in °R and T m is the mean temper- 
ature of the wall and the suspension in °R, 


Re 


s 


DVp 

s's 




(1 + v) Re g 


and 


Pr, 


f C ) 

1 P/ s M g _ 1 + 577 pr 
k g 1 + 77 g 


If equation (21) is divided by Nu (eq. (12)) and Pr is taken as 0. 75, then 

& & 

— = 16. 9 Re" 0, 3 (1 + 6ij)°- 45 (25) 
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Equation (25) is similar to Schluderberg’s equation (19) except that it also includes a 
Reynolds number effect; for example, for an average Re^ of about 25 000, the factor 
16. 9 Re“°‘ 3 =* 0. 81, and equation (25) approaches equation (19). 

The other equations are all based on using the physical properties of the suspension 
in standard pure -fluid heat- transfer correlations, although no attempt was made to in- 
clude suspension thermal conductivity and viscosity, and therefore offer no apparent ad- 
vantages. Closer examination reveals that all the equations except equation (23) can be 
reduced to the same form for h g /hg as that in equation (25). 

Gorbis and Bakhtiozin. - The heat-transfer characteristics of a gaseous suspension 
of graphite particles in vertical flow have also been studied by Gorbis and Bakhtiozin 
(ref. 6). These authors applied the Reynolds analogy to both the gas and solid phases of 
the gas-solid suspension and then superimposed the results to obtain an approximate ex- 
pression for the suspension heat-transfer coefficient. For the gas phase, the Reynolds 
analogy requires that 


_ f s 


Wg v g 


(26) 


Similarly, by definition 


Wp 



(27) 


where fg and fp are the Fanning friction factors due to the gas and particles, respec- 
tively. The appearance of p g rather than p p in the Reynolds analogy for the solid 
phase is due to Gorbis and Bakhtiozin* s definition of f p , which is based on the gas -phase 
density and is defined by the equation 


f = ap p 5 £c 

P ~ 2 % L v 2 p 


(28) 


where AP p is the additional pressure drop caused by adding particles to the gas stream. 

The summation of equations (26) and (27) with the assumption of no slip between the 
solids and the gas (V = V ) yields 

& Jr 
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L-^Vg f g 


(29) 
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so that, from equations (2) and (26), 


h /f\ 

= i + 6( JL \ 

h \f I 

g \g/ 


For the suspension by definition 
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Thus, 
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(30) 


(31a) 


(31b) 


(31c) 


(3 Id) 


The pressure drop ratio can be related to the loading ratio ?} by the simple relation 
Gasterstadt (ref. 20) proposed as far back as 1924: 


AP^ 

AP 


= i + Fn 


g 


(32) 


where F is not a constant but a complicated function of gas Reynolds number, particle 
size, etc. Combining equations (30), (31d), and (32) then gives 


h s 

= 1 + F6i] 
h„ 


(33) 


Gorbis correlated his experimental data in the form suggested by equation (33) 
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assuming F to be a function of two dimensionless groups: the gas Reynolds number 
Re^ and the particle Reynolds number Rep, where Re p = D pVp S Pg/Mg and Vp g is the 
terminal settling velocity of the particle. The terminal settling velocity is defined as 
the velocity of the particle for which the drag forces are exactly balanced by the gravita- 
tional forces and can be calculated from standard drag-coefficient correlations (e. g. , 
ref. 21). His final correlation is 

1+ ^6. 3 Reg 0 ’ 3 Re" 0, 33 j (34) 

hg \ / 


so that 


F = 


6. 3 Re 


-0.3 

g 


Re" 0 ' 33 


Equation (34) also indicates that the gas Reynolds number has less of an effect on the 
heat-transfer coefficient of a gaseous suspension than on the heat-transfer coefficient of 
a pure gas (see eqs. (15), (17), and (25)). 

The physical characteristics and size of the particles are accounted for by the di- 
mensionless group Rep, which shows that a reduction in the size of the solid particles 
will increase the heat-transfer coefficient for the same settling velocity. The correlation 
given by equation (34) applies in the range of variables 7000 < Re g < 65 000, 

5 < Re p < 800, 1 < i] < 32, 0. 84 < 6 < 1. 03, and 12 < D/D p < 133. 

In their paper, Gorbis and Bakhtiozin claim to have correlated the data of Farbar 
and Morley; however, it should be noted, that Gorbis in using his correlation assumed a 
value of 150 microns for the diameter of the particles (the smallest particle size for 
which his correlation is applicable) used in the Farbar and Morley work, whereas the 
actual average diameter of the particles was about 50 microns. Using the correct 
particle diameter (50 /jl) in his correlation does not bring close agreement with the data 
of Farbar and Morley. This lack of agreement could indicate that Gorbis’ correlation 
should not be extrapolated to particle sizes smaller than 150 microns. 

Tien . - An analytical study of heat transfer from a turbulently flowing fluid-solid 
mixture in a pipe was done by Tien (ref. 22). By making certain assumptions, Tien was 
able to solve the energy equations for both the solid and the fluid analytically. Tien’s 
analytical solution indicates that the effect of the solid particles on the heat-transfer 
coefficient is governed by the factor 6?]. This fact was also shown in the correlations 
presented by Gorbis (ref. 6) and Schluderberg, et al. (ref. 18). Tien's results are ap- 
plicable only to solid loading ratios less than one and show that in this range the heat- 
transfer coefficient varies linearly with loading ratio. The actual increase in heat- 
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transfer coefficient, however, is so slight that the advantages of adding solids in this 
range are negligible. 

Tien and Quan. - In order to verify Tien’s analytical solution (ref. 22), Tien and 
Quan (ref. 23) carried out an experimental heat-transfer study using air, and 30- and 
200-micron glass and lead particles. Gas Reynolds numbers were set at 15 000 and 
30 000, and loading ratios were varied only between 0 and 3. The experimental results 
showed a peculiar feature: the Nusselt number first decreased and then increased as the 
loading ratio was increased for fixed gas Reynolds number and particle size. This fea- 
ture was not indicated by Tien’s analytical study, and the authors proposed that the de- 
crease in Nusselt number is due to the distortion of the gas flow field by the presence of 
solid particles that was not accounted for in the theory. 

Tien and Quan did not present a correlation of their data. Qualitatively, however, 
they found that the heat-transfer coefficient ratio decreases with an increase in gas 
Reynolds number, which is in agreement with Farbar and Morley, Gorbis and Bakhtiozin, 
Farbar and Depew, and Danziger. They also found that the heat-transfer coefficient 
ratio increases with 6, as predicted by Tien’s analytical work and Gorbis’ correlation, 
and decreases with an increase in particle size, as was found by Gorbis, and Farbar and 
Depew. 

Tien and Quan also measured the maximum temperature difference between particles 
and fluids and found this to be about 20° F for the 200-micron particles, but only 3° F 
for the 30-micron particles. Thus, small particles are likely to have a negligible tem- 
perature lag, whereas a suspension of large particles is likely to be far from equilibrium 
conditions. 

Mickley and Trilling. - Mickley and Trilling (ref. 24) conducted their research to 
determine the effect of the presence of solids on heat-transfer conditions at surfaces in 
contact with a fluidized gas -solid mixture. The experimental work was carried out in 
vertical fluidized beds 1 to 4 inches in diameter. The fluidized mixture consisted of 
glass spheres with particle diameters ranging from 40 to 450 microns suspended in an 
upward flowing stream of air at superficial velocities varying between 0. 8 and 15 feet 
per second. Several groups of particles of fixed size were tested in this range. 

The work of Mickley and Trilling was carried out in fluidized beds, and hence the 
direct validity or quantitative applicability of their results is questionable for use with 
suspensions; however, the important trends and variations observed in their data might 
be qualitatively applicable to solid suspensions flowing in tubes. 

Mickley and Trilling found that for the four groups of large-size particles (150 to 
450p), the slope of the line that resulted when the log of the heat-transfer coefficient 
was plotted against the log of the particle concentration was relatively constant at a value 
between 0. 46 and 0. 55. This value is in good agreement with the suspension work of 
Farbar and Morley, Danziger, and Schluderberg. They also found that at any given 
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concentration, the heat-transfer coefficient varies at a rate inversely proportional to the 
0. 61 to 0. 74 power of the particle diameter for the larger sized particles. However, 
this effect of particle diameter on heat-transfer coefficient decreased considerably for 
the smaller sized particles (D p < 150 n). This behavior was attributed to agglomeration 
of the smaller sized particles. 

List . - The final heat-transfer study that will be reported in this technical note was 
done by List (ref. 25). List determined heat-transfer coefficients for four sizes of glass 
beads and two sizes of copper shot transported by air in a vertical heating section. The 
mixtures were electrically heated with an internal rod heating element. The solid flow 
rate for a given run was set by choosing one of a set of cone-shaped nozzles, precali- 
brated by determining the average weight of the solids delivered under static conditions. 
The smallest bead used was 51 microns in diameter, about the same size as the catalyst 
used by Farbar and Morley and by Danziger. List's data showed that h g was inversely 
proportional to the solid particle size to the 0. 04 to 0. 09 power. At low velocities, the 
effect of particle size on h g was found to increase. Mickley and Trilling found that h g 
varied inversely proportional to the 0. 61 to 0. 74 power of particle size. Indications are 
therefore that as the velocity increases from that of a fluidized bed to that of suspension 
flow, the effect of particle size becomes smaller. 

Although List reports extensive data that would be valuable to any heat-transfer 
study, Danziger (ref. 5) believed List's data to be unreliable because the solid inlet 
nozzles were calibrated at static conditions rather than at dynamic conditions. There- 
fore, the data of List have not been used in the attempt to find an overall correlation. 


Correlations 

From the literature cited, it is apparent that the important parameters that deter- 
mine heat transfer in a gas-solid suspension are the loading ratio r\, the ratio of speci- 
fic heat of the particles to the specific heat of the gas 6, the particle diameter Dp, and 
the gas Reynolds number Re g . The effect of each of these parameters on the h g /h g 
ratio will now be studied in more detail based on the results of the aforementioned 
studies, and an attempt will then be made to establish a consistent workable correlation. 

Effect of loading ratio 7? and specific heat ratio 5. - The heat-transfer correlations 
deduced in the references cited are summarized below in terms of the heat-transfer 
coefficient ratio: 

Farbar and Morley: 


K ^.= 6.8 Re -0 ' 2 ?; 0 ’ 45 
g 


(15) 
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Danziger: 


i=3.7Re; 01 V- 45 

\ g 

Schluderberg: 

^.= 0.18 (1+ Oi?)°* 45 

h 


(17) 


(19) 


The Franklin Institute: 


^ = 16. 9 Re"°* 3 (1 + &n)°' 45 

Gorbis and Bakhtiozin: 

— = 1 + ((5. 3 Re" 0 * 3 Re" 0 * 33 ) &n 

K V S p / 


(25) 


(34) 


The one fact that finds unanimous agreement from these equations is that the heat- 
transfer coefficient ratio increases as the loading ratio increases although three different 
forms of the loading ratio dependence are suggested. If all of the other parameters of the 
system are held constant, these forms may be written as 


A A __ 

— = k-,- a 




(35) 
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— = k 2 (l + kg ??) 1 
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(36) 


h s 

— = 1 + k 4 ?? 

h„ * 
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(37) 


where the k's are, in general, functions of other variables, and a and b are positive 
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exponents less than one. 

In the range of loading ratios of primary interest (77 < 10), there has been little ex- 
perimental work done to determine the effect of 6 on h g /hg although it almost certainly 
exerts a strong influence on the heat-transfer coefficient ratio. The only work where 6 
has been specifically investigated experimentally is that of Tien and Quan at low loading 
ratios. Their work qualitatively indicates that hVh increases with an increase in 6. 
Tien’s analytical work indicates that the factor 5rj should be treated as a unit, which 
also is predicted by the Reynolds analogy as presented by Gorbis. Schluderberg and the 
Franklin Institute also use die factor 677 as a unit in their correlations. Although the 
available data are very limited as to the effect of 6, it seems reasonable to assume that 
the factor 677, which can be thought of as a weighted specific heat ratio, should appear as 
a unit. Thus, equations (15) and (17) are rewritten for the respective gases and particles 
used: 

— = 6. 4 Re' 0 ' 2 («rj) 0 ' 45 (38) 

h g 

and 

— » 4. 0 Re -0, 14 (5t7) 0, 45 (39) 

h 

g 

\ 

It should be noted at this point that some disagreement exists concerning the value of 
6 of the silica -alumina catalyst used by both Danziger and by Farbar and Morley 
(table I). The procedure utilized in the preceding paragraphs, however, negates any 

possible errors introduced by an incor- 
rect 6 in obtaining their respective 
constants since the constants appearing 
in equations (38) and (39) are now pre- 
sumably independent of 6. 

In order to see how closely the cor- 
relations agree with one another, the 
heat-transfer coefficient ratio as pre- 
dicted by equations (38), (39), (25), and 
(34) was plotted against 677 in figure 6. 
The values are plotted for a gas Reynolds 
number arbitrarily set at 20 000 and for 
the smallest allowable particle diameter 



Weighted specific heat ratio, (nj 

Figure 6. - Effect of weighted specific heat ratio on heat-transfer 
coefficient ratio. Gas Reynolds number, 20 000. 
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Figure 7. - Effect of weighted specific heat 
ratio on heat-transfer coefficient ratio. 
Gas Reynolds number, 20 000. 


(150 fi) in the Gorbis and Bakhtiozin correlation. The 
figure shows close agreement between all the correl- 
ations indicating that the form (eqs. (38) and (39)) 

— = kj (&?) 0 ' 45 (40) 

h o 

is a good representative correlation. On the other 
hand, a plot of (h g /hg) - 1 against 617 , obtained from 
the correlation equations as presented in figure 7, 
also shows fairly good agreement between all the cor- 
relations. This agreement indicates that an equation 
of the form 


h 

— = 1 + k’ 5f] 
h 4 

g 


(41) 


is also a good representative correlation in the same range. Thus it seems that either 
equation (40) or (41) is applicable in the range of interest 2 < 77 < 10. However, while 
h s /hg in equation (41) reduces to unity as it should when 77 — 0, equation (40) does not. 
Equation (41), however, gives high values for h /h when compared to Danziger's ex- 

S 5 ' 

perimental data at loading ratios greater than 15, whereas equation (40) is in agreement 
with Danziger's data in the entire loading ratio range. 

The form of equation (25) 


— = k ’ 2 (1 + 677 )°’ 45 (42) 

h g 

seems to offer no advantage over the other two forms and will be abandoned. 

Effect of particle diameter . - Another parameter of importance on which little ex- 
perimental or analytical work is available is the particle diameter. It is generally ac- 
knowledged that there should be some decrease in the heat-transfer coefficient ratio as 
the particle diameter is increased. The experimental works that attempted to find the 
particle diameter dependence are the work of Gorbis and Bakhtiozin, Farbar and Depew, 
List, and the work of Mickley and Trilling with fluidized beds. 

Gorbis and Bakhtiozin expressed the particle diameter dependence in terms of the 
dimensionless group Re p (eq. (34)). This group can be expressed specifically in terms 
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Figure 9. - Effect of particle size on ratio of heat-transfer coefficients. Gas Reynolds number, 
26 500; weighted specific heat ratio, & 


of particle diameter by eliminating Vp g explicitly with the help of standard drag coef- 
ficient data (ref. 21) and shows that particle diameter exerts a small influence on h g /hg. 
List's work (ref. 25) indicates that the heat-transfer coefficient ratio varies inversely 
with the 0. 04 to 0. 09 power of the particle diameter. Farbar and Depew show a much 
larger effect of particle diameter and indicate that with particles of 200 -micron diameter 
the heat -transfer coefficient of a suspension is essentially the same as that of a pure gas. 
Mickley and Trilling (fluidized beds) suggest that h g /hg varies inversely with the 0. 6 to 
0. 7 power of the particle diameter for particles greater than 150 microns and that the 
effect is smaller for smaller particles. This contention is in contradiction with the re- 
sults of Gorbis and Bakhtiozin and List; however, this large particle diameter depend- 
ence is believed attributable to the low velocities encountered in fluidized beds. 

A plot of h s /hg against D p at a constant gas Reynolds number of 26 500 is pre- 
sented in figure 8 for a value of 5t] equal to 2 and in figure 9 for a value of 6?j equal 
to 8. The correlations of Farbar and Morley, Danziger, and Schluderberg (or the 
Franklin Institute) appear as single points on these plots since these correlations were 
each based on only one particle diameter. List's work is represented by a line of aver- 
age slope -0.065 and was arbitrarily passed through Danziger's result since List's 
smallest particle size was 50 microns and he presented no general correlation for his 
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work. Far bar and Depew's experimental results and Gorbis's correlation are also pre- 
sented in the figures within their ranges of applicability. It should be pointed out that 
Gorbis's correlation, in particular, should not be extrapolated to particle diameters 
below 150 microns especially at high 6/j since the correlation shows a rapid rise in 
slope at this point, which is inconsistent with the experimental data in the low particle 
diameter range. As can be seen from the figures, with the exception of Farbar and 
Depew’s data, at high &q the effect of particle diameter is indicated to be small in the 
entire particle diameter range from 5 to 1000 microns. 

For application in a Brayton cycle and other pumped circuits the particle diameter 
should be kept as small as possible to prevent surface erosion and any appreciable lag 
between the particle and gas in the rotating equipment; yet the particle diameter must be 
large enough to avoid stability problems. Stability problems arise if the particles sepa- 
rate from the gas stream and adhere to the boundary surface of the heat-transfer equip- 
ment because of a thermal gradient. For example, a visual inspection of Schluder berg's 
experimental loop showed that the gas suspension was remarkably free of erosion, plug- 
ging, or particle adherence to system surfaces except on the cooling surfaces of the loop 
where graphite particles were found to plate out. The instability of the suspension in 
these areas was caused by an anisotropic Brownian particle motion in the direction of 
the decreasing temperature gradient. Recent information obtained from the Bureau of 
Mines (ref. 26) indicates that the use of spherical graphite particles larger than 10 mi- 
crons will minimize the stability problem. Thus, the range of particle diameters of in- 
terest for Brayton cycle and related work is of the order of 10 to 50 microns. In this 
range there appears to be no appreciable effect of particle diameter on the ratio h g /hg- 

Effect of gas Reynolds number. - Most of the literature indicates that a decrease of 
hg/hg should result when there is an increase in gas Reynolds number. This decrease 
can be explained by a reduction in the relative turbulence -producing effect of the motion 
of the solid particles in the gaseous flow, the turbulence of which increases with an in- 
crease in the Reynolds number. 

Plots of h g /h against gas Reynolds number for two representative values of &q 
of 2 and 10 are shown in figures 10 and 11, respectively. The figures show close agree- 
ments between the correlations of Farbar and Morley (eq. (38)), Danziger (eq. (39)), 
Schluderberg (eq. (25)), and Gorbis and Bakhtiozin (eq. (34)) for 150 -micron particles, 
and the data of Farbar and Depew for 30-micron particles. The data of Farbar and 
Depew does not extend to dq = 10 and therefore could not be plotted in figure 11. A 
representative average slope of these curves can be taken as -0.21, which is in good 
agreement with that predicted by Farbar and Morley' s correlation, as shown by the 
faired dashed lines in figures 10 and 11. 

Good agreement among the various correlations can also be obtained in terms of 
(h g /h ) - 1, as shown by figure 12. Figure 12 is a plot of (h/h ) - 1 against Re for 

° 5 ° 5 o 
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Figure 10. - Effect of gas Reynolds number on heat-transfer coefficient ratio. Weighted specific heat ratio, 2. 


dri = 10. An average representative slope of these curves, as indicated by the faired 
dashed line, is -0. 32, which is in good agreement with Gorbis's correlation. 

Final form of correlation. - The final forms of equations (40) and (41) based on the 
representative (dashed) lines drawn in figures 10, 11, and 12 appear as follows 


-1 = 7.6 Re: 0,21 (677) 0,45 

h & 


(43) 


-1=1 + 4. 0 Re -0 ' 32 577 


(44) 


These equations have been plotted against in figure 13 for three different gas 
Reynolds numbers. As can be seen from the figure, the agreement between the correla- 
tions in the range of 2 < t? < 10 is good. 
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Figure 12. - Effect of gas Reynolds number on heat-transfer coefficient ratio factor. Weighted specific heat ratio, 10. 
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Figure 13. - Comparison of two final heat-transfer coeffi- 
cient correlations at various gas Reynolds numbers. 


flow in a pipe is given by the Fanning equation 


These correlations should be applica- 
ble in the range of gas Reynolds numbers 
from 5000 to 100 000, 6 rj between 2 and 10, 
and to particle diameters between 5 and 
50 microns. As indicated in the preceding 
section, the particle diameter effect in 
this range has been regarded as negligible. 

PRESSURE DROP 

hi view of the widespread use of gas- 
solid conveying in a variety of industries, 
it would be expected that the friction and 
pressure drop flow phenomena would be 
well defined and understood. Investigation 
of the literature, however, indicated that 
the prediction of pressure loss in gas -solid 
suspensions is still an empirical art. In 
the ensuing sections an attempt will be 
made to find a working correlation that will 
predict pressure drop of suspension flow. 

Pure Gas Relations 

For a pure gas, the pressure drop for 


AP = 2f — p — §L 

g g D p gg c 


(31c) 


where the value of the gas friction factor f g is well established for all flow regimes. 

For fully developed turbulent flow in smooth tubes under adiabatic conditions 

f g = 0. 046 Reg 0, 2 for 10 000 < Re g < 100 000 ( 45 ) 

It would be convenient to develop an equation similar in form to the Fanning equation 
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with appropriate modifications to predict the pressure drop of a gas -solid suspension. 

The friction factor for flow of a gas-solid suspension, however, is not expected to be 
merely a simple function of gas Reynolds number but may depend on other variables such 
as loading ratio t] and particle diameter D . 

Flow Regimes 

‘ There have been a number of studies made to determine the effect of the suspension 
flow patterns (gas and solid velocity profiles) on the pressure drop of suspensions 
(refs. 27, 28, and 29). The most recent effort to determine the mechanism of gas-solid 
flow in tubes is the work of Peskin and Dwyer (ref. 30). From their experimental work 
they have concluded that gas -solid flow can be divided into the four following distinct flow 
regimes: 

(1) Particles cause viscous disturbances and increase the size of the laminar sub- 
layer and thus a decrease in the shearing stress at the wall. In this case, the particles 
are very close together but occupy a small volume. 

(2) Particles cause mostly viscous disturbances, but are too far apart to affect the 
gas velocity profile and hence the laminar sublayer. The shearing stress at the wall is 
unchanged. 

(3) Particles cause inertial disturbances that alter the gas velocity profile and de- 
creases the size of the laminar sublayer. In this case, the wall shearing stress in- 
creases. 

(4) Particles occupy a large volume and change the geometry of the flow. This situ- 
ation is similar to the flow in packed beds. 

It appears therefore that the variation of the friction factor of a gas -solid suspension 
with loading ratio and other parameters may depend on which flow regime is applicable 
and may not be generalized for all gas -solid flows. 


Analytical Models 

Several analytical models for describing the pressure loss of gas -solid suspensions 
are available based on the assumption that the mixture flows at a velocity greater than 
either the saltation velocity in horizontal tubes (the lowest possible velocity required to 
keep a specific particle in suspension) or the choking velocity in vertical tubes (the mini- 
mum velocity required to prevent the suspension from being transported up the tube in 
slug flow). Under these conditions the pressure drop for the general flow of gas -solid 
suspensions through tubes is considered to be composed of the following six forces: 
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(1) The friction of the gas against the pipe wall 

(2) The friction between the solids and the pipe wall 

(3) The drag force required to move the solids through the pipe 

(4) The force required to accelerate the gas to its equilibrium velocity profile 

(5) The force required to accelerate the particles to their equilibrium velocity 

(6) In vertical tubes, the force required to support the weight of solids and gas 
The following models were employed to determine pressure drop in a gas -solid suspen- 
sion and are based on the inclusion of the preceding forces in varying degrees and forms. 

Drag coefficient model . - Zen z and Othmer (ref. 21) present the following equation ' 
proposed by Hinkle (ref. 31) for predicting the pressure drop of a dilute suspension: 


4P S = 


2g, 
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(46) 


Fluid 


Particle Fluid to Particle 


accelera- accelera- pipe 
tion tion friction 


friction 

and 

drag 


Weight of 
particles 
(Vertical 
tubes only) 


where 


p* density of solid, p - p , lb solid/ff 5 mixture 

r ° 5 


f* particle friction factor based on particle velocity and particle density in mixture 

r 


For a horizontal tube of constant diameter that is of sufficient length so that acceleration 
forces can be neglected (particles and fluid having reached their respective equilibrium 
velocities), this equation reduces to 
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or with the use of equation (31c) 
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But since t] is the ratio of the mass flow rate of particles to the mass flow rate of pure 
gas 


W„ p'V 

tj = — = 1E_P 

W_ p v 
g g 


then 
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f V 
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V 


(49) 


(50) 


Equation (50) is simply another form of the Gasterstadt relation (eq. (32)) where 
F = fpVp/fgVg. Before equation (50) can be used to determine AP g , however, a means 
of obtaining both fp and Vp is necessary. 

For particle sizes showing a significant velocity lag or slip between the particles 
and the transporting fluid, Hinkle (ref. 31) assumed the pressure drop due to the pres- 
ence of the particles to be caused entirely by the drag of the fluid on the particles. In 
this case, the particle -to -wall friction can be neglected in comparison with the drag. 
Based on this assumption, the pressure drop due to the particle drag is 


2 *>; v p l . nc p p B (v g - v 2 ,d p l 

g c D 8g c 


or 


, _ nC D Pg 
P 16 pj, 



,d 2 d 


(51) 


where n is the number of particles per cubic foot of mixture and Cjj is the standard 
drag coefficient for spherical particles defined by 


C D = 


8F P g c 

^ g (v g - V 2 


where F D is the drag force and - Vp is the relative velocity between the gas and 
particle. 
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Since 


n 
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then 
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P p p p 


(52) 


In order to employ this model the particle velocity must be known. If experimental 
values of V p are not available, the particle velocity can be crudely estimated by assum- 
ing V to be equal to V - V „ where the terminal settling velocity V„ 0 can be ob- 
p g ps ps> 

tained from the standard drag coefficient - Reynolds number correlation (ref. 21): how- 
ever, the attendant difficulty of obtaining an accurate particle velocity makes this model 
difficult to apply in obtaining the particle friction factor. In addition, the effect of 
particle-to-wall friction has been entirely omitted. For example, for particles small 
enough to show no slip (Vg = V p ), this model would predict no additional pressure drop 
due to the particles. 

Equivalent friction factor model . - The inherent disadvantages of the previous model 
indicate the desirability of having a simpler model to predict suspension pressure drops. 
The equivalent friction factor model simply defines the pressure drop of the suspension 
as 


V 2 

AP s = 2f W - — 
S g D g c 


(31a) 


where f^ is an equivalent suspension friction factor and is the density of the gas in 
the suspension. If this equation is divided by the Fanning equation for the same gas con- 
ditions in the tube then 


AP 

AP 



(3 Id) 


An alternative definition of the equivalent friction factor based on suspension density 
is given by 


32 



( 53 ) 


V 2 

AP s = 2i s p s ~ ~ 
S S ® D g c 


Thus f ' is related to f by the relation 
s s 


's = *s~ = f s (I + ’» 
Pa 


(54) 


In both cases, f_/f_ as well as f' /f„ are presumably functions of loading ratio, gas 
Reynolds number, particle diameter, and other properties of the suspension that must be 
determined experimentally. 

Reynolds analogy. - Another approach to the determination of the suspension friction 
factor is the use of the Reynolds analogy. The Reynolds analogy is of historical impor- 
tance as the first recognition of analogous behavior of momentum- and heat-transfer 
rates. The Reynolds analogy can be applied to gas-particle suspensions if the suspension 
is assumed to behave as a homogeneous fluid and if an accurate relation for the heat- 
transfer coefficient of suspensions is assumed available. 

By applying the Reynolds analogy to the gas and solid phases separately, as was done 
by Gorbis in reference 6, the relation between f p /fg and h g /hg is obtained: 



1 + 


£ g 


(30) 


By definition in this development 


f = f + f 
s g p 


(31b) 


where f' is defined by equation (31a), f is defined by equation (28) , and f is defined 
s p & 

by equation (31c). Thus with the use of equation (30) 



s 

g 



(55) 


The Reynolds analogy can also be applied by comparing the gas directly with the 
suspension rather than by superposition of the gas and solid phases as was done pre- 
viously. The Reynolds analogy yields 
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for the gas 


(26) 


g 


- f g 




and 


h f 

2 = — for the suspension 


where f is the equivalent friction factor defined by equation (53). Dividing f by f 
s s & 

and assuming equal volumetric rates of flow so that Vg = V g give 



p. 

( C P1 

k 

f g h g 

p s l 

( c p) 

1 

s_ 


In terms of loading ratio rj for equal velocities (eq. (lc)) and ratio of specific heats 6 
(eq. (3)), 


f s _ h s 1 I + 77 
f g h g 1 + rj 1 + 6»] 


or 


t„ h_ I +677 
g g 

With the use of equation (54), 

= jb l_+n 

f h_ I +677 
g g ' 


(56) 


(57) 


Thus, either equation (55) or (57) can be used to relate the equivalent friction factor of 
the suspension to the heat-transfer coefficient of the suspension, depending on the method 
of application of the Reynolds analogy. When 5=1, both forms reduce to the simple 
expression 
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( 58 ) 



Thus, once a workable relation for h g /hg is accepted, a reasonable relation for f g /f 
might also result. 

Eddy viscosity model . - In a recent paper, Julian and Dukler (ref. 32) attempted to 
find a correlation for pressure drop caused by a gas -solid suspension by use of an eddy 
viscosity model. In their work, they suggest that for dilute -phase transport the solids 
make their presence felt primarily by modifying the local turbulence in the gas phase, 
increasing the turbulent fluctuations, mixing length, and eddy viscosity, and consequently, 
frictional pressure drop. 

Julian and Dukler's analysis is based on a modification of Gill and Scher's expres- 
sion (ref. 33) for the eddy viscosity e* of a pure fluid flowing through a tube. Gill and 
Scher's equation is 


e' 


= K 2 y 2 



(59) 


where K is the universal von Karman constant for pure gas originally reported as 0. 4 
and found by Deissler (ref. 34) to be closer to 0. 36. Julian and Dukler modified the Gill 
and Scher equation to take into account the effect of particles by redefining. 

K = k(l + r]) m (60) 

The unknown constants k and m are to be determined from experimental data. By a 
trial-and-error procedure, using the experimental data of six different investigators 
(refs. 31 and 35 to 39) Julian and Dukler were able to show that a log plot of K against 
l + ?7 is essentially a straight line as predicted by their model. For loading ratios below 
12 they found that 


K = 0. 36(1 + 77 )°* 25 


(61) 


and for loading ratios above 12 


K = 0. 11(1 + 7j)°" 8 (62) 

Equation (61) is especially encouraging since at tj = 0, K reduces to 0. 36, which is the 
accepted value for flow of a particle -free gas. 
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Since the eddy viscosity can be readily related to the equivalent friction factor f ' 

s 

of equation (31a) (ref. 32), Julian and Dukler’s model can be used to obtain suspension 
pressure drops. These results are discussed in a later section (Correlations section). 

Data Sources 

In order to evaluate fully the validity of the aforementioned models and mechanisms, 
an extensive literature search was undertaken. The purpose of this literature search 
was to uncover experimentally obtained data and correlations capable of being used in 
this evaluation procedure. A brief review of the literature available in this field is pre- 
sented in the following sections, and a summary of the operating conditions and the phys- 
ical properties of gases and solids used in each experiment is presented in table n. The 

TABLE n. - SUMMARY OF OPERATING CONDITIONS USED IN PRESSURE DROP WORK CITED 


Particle 

Gas 

Particle 

diameter, 

D P’ 

Loading 

ratio, 

V 

Gas 

velocity, 

> 

ft/sec 

Equivalent 

tube 

diameter, 

D, 

in. 

Gas Reynolds 
number, 

Re g 

Tube 

length, 

L, 

in. 

Source 

Glass 

Air 

00 

p 

>-» 

I-* 

o 

0 to 4 

100 

3 

132 000 

36 

Peskin and 
Dwyer 
(ref. 30) 

Glass 

Air 

80, 110 

0 to 14 

60 to 100 

3 

100 000 to 
150 000 

36 

Peskin 
(ref. 40) 

Glass 

Air 

36, 97 

0. 7 to 7 

10 to 90 

0.5 

10 000 to 
60 000 

158 

Mehta, Smith, 
and Comings 
(ref. 41) 

Alumina - 

silica 

catalyst 

Air 

10 to 220 
(avg 50) 

0 to 16 

50 to 150 

0.69 

17 000 to 
50 000 

87 

Farbar 
(ref. 42) 

Millet, 
peas, turnip 
seeds, pine 
kernels and 
sunflower 
seeds 

Air 

1280, 5760, 
1900, 5850, 
4000, 8400 

0 to 15 

65 to 130 

5 

110 000 to 
321 000 

80 

Dogin and 
Lebedev 
(ref. 43) 

Sand, clo- 
ver seed, 
wheat 

Air 

200, 730, 
330, 1150, 
440, 4000 

0.6 to 41. 6 

35 to 150 

0.5 

9000 to 
39 000 

300 

Vogt and 
White 
(ref. 44) 
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experiments were conducted at room temperature with pressure levels ranging from 1 to 
3 atmospheres among the various data sources. 

Peskin and Dwyer . - In order to verify the drag coefficient model (eqs. (50) and (52)), 
Peskin and Dwyer (ref. 30) measured pressure drop and both fluid and particle velocity 
profiles for loading ratios up to 4 by using 80- and 110-micron-diameter glass particles 
in air in a 3-inch-square duct. For these small, relatively widely spaced particles 
(regime 2), Peskin found that the flow disturbances caused by the particles are entirely 
viscous since the particles because of their small size do not possess enough inertia in 
themselves to cause an inertial disturbance in the gas. Also, because of the relatively 
large spacing between the particles, the gas velocity profile was found to revert back to 
a universal profile between the particles. The average velocity of the particles, however, 
showed a large deviation from the velocity calculated by using the standard drag coeffi- 
cient for the settling velocity; that is, on the basis of the experimentally measured par- 
ticle velocity, the drag coefficient calculated from the pressure drop measurements 
was much higher than that predicted from the standard correlation (eq. (52)). Peskin 
attributes the high drag coefficients to the added drag of the longitudinal component of the 
turbulence of the fluid rather than to shortcomings in the assumed model. Because of 
the disparity between the model and Peskin’ s data, the drag coefficient model does not 
appear useful for predicting suspension pressure drops. 

Peskin . - In another investigation, Peskin (ref. 40) also correlated some of his 
pressure drop data in terms of the equivalent friction factor f defined by equation (53). 
Sets of curves of f g against gas Reynolds number were obtained for two different glass 
particle sizes (80- and 110-p. diam.) at Reynolds numbers of 100 000 to 150 000 and load- 


ing ratios up to 14 in a 3-inch-square duct. The faired variations of these curves after 
conversion to f' are presented in figure 14 and show that in all cases the equivalent 



Loading ratio, 17 


friction factor f g increases with increasing 
solids loading ratio and with decreasing gas 
Reynolds number. In addition, the curves 
indicate that at a given loading ratio and 
Reynolds number, the friction factor is 
smaller for the 80-micron particles than for 
the 110-micron particles. These results are 
contrary to the heat-transfer work discussed 
earlier (p. 23), where it was shown that 
heat-transfer coefficients tend to decrease 
slightly with increasing particle diameter. 
The effects of particle diameter on heat- 


Figure 14. - Fanning friction factor of suspension as func- 
tion of loading ratio as computed from faired data of ref- 
erence 40. 


transfer coefficients, however, may be dif- 
ferent from the effects on friction factor 
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because of the velocity lags inherent with larger particles. 

It should also be noted at this point that the results for the friction factor of the 
110-micron particles are somewhat dubious in view of certain experimental inaccuracies 
in the runs for this particle size (private communication from R. L. Peskin). It is 
therefore conceivable that the large -particle -diameter effect reported by Peskin may not 
be correct. 

Mehta, Smith, and C omings. - In the work of Mehta, Smith, and Comings (ref. 41), 
36- and 97-micron-diameter glass beads were each suspended in air flowing through a 
1/2-inch-diameter pipe. The gas Reynolds number was varied between 10 000 and 
60 000, and the loading ratio was varied between 0. 7 and 4. 7 for the 36 -micron particles 
and between 0. 7 and 7. 0 for the 97-micron particles. As a result of their study, the 
authors concluded that the pressure drop in air-solid transport systems is dependent 
upon the type of particle flow. The authors of reference 41 postulated that the 97-micron 
solids were primarily in ’’bouncing'’ flow, where the particles are in unsteady motion 
and frequently collide with the wall of the pipe. The 36-micron particles were assumed 
to be in suspension flow, where the particles flow in a suspended condition that is main- 
tained by the finite slip velocity between the particles and the gas. 

The authors attempted to correlate their data using the drag coefficient model 

(eq. (52)) with no success. They finally proposed a correlation based on a so-called 

mixture friction factor f defined by 

m 


AP 


s 


f — V?p 


2g r 



(63) 


where a’ is a constant chosen by the authors as 0. 3 for the 36-micron particles and 
1. 0 for the 97-micron particles. Upon plotting f from their pressure drop data 
against gas Reynolds number, they found f to be roughly constant, or independent of 
solid flow rates, for both particle sizes at Reynolds numbers greater than 30 000. The 
average values of f were 0. 016 for the 36-micron particles and 0. 035 for 97-micron 
particles. At gas Reynolds numbers less than 30 000 these values increased as Re^ 
was decreased. 

In Mehta's experiments, V ^/Vg f° r the 36 -micron particles remained relatively 
constant at about 0. 7, whereas for the 97-micron particles, V^/Vg varied considerably 
(from about 0. 22 to 0. 67). The reason for this might be that the 97-micron particles 
were not fully accelerated when the pressure drop measurements were recorded, and 
this could therefore account for the much higher pressure drop observed with the 
97-micron particles. 
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Dividing equation (63) by equation (31a) gives the following relation between fg and 


m 



(64) 


Substituting the constant values for f m suggested by the authors for each particle size 
gives 


f 


T 

s 


= 0. 004 + 0. 004 



V 



for 36 -micron particles 


(65) 


and 


f’ = 0. 0088 + 0. 0088 -E ij for 97-micron particles (66) 

S V„ 

g 


The correlations given by equations (65) and (66) are difficult to apply, however, inas- 
much as they depend on the magnitude of the particle velocity, which is a difficult quan- 
tity to determine. In addition, the correlation is only applicable for Re^ > 30 000 since 
at lower gas Reynolds numbers, f was not found by the authors to remain constant. 



Figure 15. - Fanning friction factor of suspension as func- 
tion of loading ratio as computed from data of reference 41. 


A plot of f_ against t] for several 
values of Reg calculated from the Mehta, 
Smith, and Comings original horizontal- 
flow data (for 77 > 1) is presented in fig- 
ure 15. The curve deduced from their 
36-micron particle correlation according 
to equation (65), with V p/V g = 0. 7 and 
Reg greater than 30 000, is included as a 
comparison (denoted by the dashed line). 
Faired lines through the data points indi- 
cate that f’ varies with approximately the 
0.4 to 0. 5 power of 77 , whereas the cor- 
relation gives a smaller dependence on 77 . 

The figure also shows that f' decreases 

s 

with increasing gas Reynolds number. 

These slopes are in general agreement with 
the data of Peskin shown in figure 14. 
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friction factor of suspension to friction factor of gas (data 
from ref. 42. ) 


Farbar . - Farbar (ref. 42) measured 
the flow characteristics of solid-gas mix- 
tures in both a horizontal and a vertical 
pipe 0. 69 inch in diameter using a mixture 
consisting of silica -alumina catalyst and 
air. The loading ratio was varied between 
0 and 16, and the particles had a size dis- 
tribution varying from less than 10 microns 
to greater than 220 microns with an average 
diameter of 50 microns. The gas Reynolds 
number was varied between 17 000 and 


50 000. 

As a result of his studies, Farbar was able to report that AP g / AP^ = f g /f^ ap- 
peared to decrease when the gas flow rate Re was increased. Farbar also reported 

that AP/AP increased as 77 was increased up to 77 = 10. Although his data above 
s g 

7] = 10 were inadequate to establish a definite variation, he proposed that above this 
value the loading ratio no longer has an effect on the pressure drop of the suspension. 

It should be noted, however, that Farbar presents no correlation for his data and that the 
data (as given in his paper) for vertical transport are inconveniently presented because 
the author neglected to subtract the gravity head (the pressure drop caused by support- 
ing a vertical column of suspension) from the total pressure drop to obtain the pressure 
drop caused by friction alone. 

A plot of the curve given in Farbar’ s paper to represent the horizontal tube pressure 
drop data appears in figure 16. Although Farbar chose to draw a single line through all 
of his data, several lines could have been drawn for different values of Re . The curve 

T g 

shows, however, that f^/f varies approximately with the 0. 42 power of 77 which is in 
good agreement with both the data of Mehta and Peskin. 

Dogin and Lebedev. - In the work of Dogin and Lebedev (ref. 43), an attempt was 
made to determine the dependence of the friction factor f’ of gas-solid suspensions on 
the loading ratio, flow velocity, specific weight, and dimensions of the particles being 
conveyed. The materials conveyed by air in a 5-inch-diameter horizontal pipe were 
millet (Dp = 1900ju), peas (D p = 5760 11 ), turnip seeds (D p = 1280ju), wheat (D p = 4000 /i), 
pine kernels (D^ = 8400 11 ), and sunflower seeds (D p = 5850 n). The loading ratio was 
varied between 0 and 15, and the gas Reynolds numbers ranged from 110 000 to 321 000. 
Dogin and Lebedev correlated their data by using the form presented by Gasterstadt 


— = 1 + Ft] (32) 

f g 
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They observed that fg/fg decreased as the flow velocity of the gas increased and found 
that (fg/fg) - 1 was proportional to the ratio of the particle diameter to the pipe diameter 
taken to the 0. 1 power and to the ratio of the densities of particle and gas Pp/Pg- The 
results of their investigation were summarized in their final correlation, which was 
given as 




0 . 1 0 . 4 

(Re)^ p 

-H 7] 


Fr 


0.5 


(67a) 


9 

where Fr is the gas Froude number equal to Vg/gD, and A is a parameter that 

varied, depending on pipe roughness, between 1X10“° and 2.2x10“ in the experiments of 
Dogin and Lebedev. 

Substituting equation (45) for f in equation (67a) indicates that for a given gas, 

g 

particle density, and tube diameter 


f T 

-= 1 + C(D p )°- ^Vg)' 0 ’ \ (67b) 

f g 

where C is a constant depending on the system. Equation (67b) shows that (fl/f_) - 1 

s g 

should increase slightly with particle diameter and increase with decreasing gas veloc- 
ity. These trends are in good agreement with most of the literature. 

Vogt and White . -In a paper by Vogt and White (ref. 44), the pressure differential 
required to produce steady flow of suspensions of sand (D p = 200, 330, 440, and 730p), 
clover seed (D p = 1150ju), and wheat (D p = 4000;ii) in air through vertical and horizontal 
1/2-inch -diameter pipes is presented. In their work, the loading ratio was varied be- 
tween 0. 6 and 41. 6, and the gas Reynolds number ranged from 9000 to 39 000. 

Vogt and White correlated their data in a form similar to that of Dogin and Lebedev 




( 68 ) 


where B and f are functions of the dimensionless group Re (Cjj)*'^. The large 
dependence of f’ /f., on D/D in equation (68) has been disputed by Belden and Kassel 
(ref. 45) and others. Clark, et al. (ref. 38) reported that the sand rims might have been 
influenced by high electrostatic charges, which could have been responsible for the high 
pressure drops reported for these particles. Vogt and White's correlation and also all 
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TABLE m. - RANGE OF DATA USED BY JULIAN (REF. 32) FOR EVALUATION OF HIS MODEL 


Particle 

Gas 

Particle 

diameter, 

d p’ 

Loading 

ratio, 

V 

Tenite, 

polystyrene, 

aluminum 

Air 

2290 to 
6360 

0. 5 to 4. 5 

Glass 

Air 

25 to 200 

0.5 to 3.5 

Cress 

seed 

Air 

1000 

2.2 to 7.4 
0.4 to 12.1 

Sand, glass, 
polystyrene 

Air 

435 to 
940 

Sand, crack- 
ing catalyst 

Air, 

co 2 

200 to 
510 

0.5 to 36 

j 

Glass, 
fertilizer, 
cress, brass, 
aluminum 

Air 

125 to 
1525 

1.6 to 16.5 


Tube 

inside 

diameter, 

D, 

in. 

Gas Reynolds 
number, 

Re g 

Source 

2,3 

178 000 

Hinkle 
(ref. 31) 

0.7 

13 500 to 
27 400 

Depew 
(ref. 39) 

1 

39 000 to 
75 000 

Clark, et al. 
(ref. 38) 

1.5, 2 

33 000 to 
59 000 

Helander 
(ref. 35) 

0.25, 

3700 to 

Hariu and 

0.5 

11 000 

Molstad 
(ref. 37) 

1 

43 000 to 
78 000 

Mitlin 
(ref. 36) 


of his runs with sand are therefore highly questionable. 

Other work . - Other important experimental studies involving pressure drop of gas- 
solid suspensions include the work of Hariu and Molstad (ref. 37); Clark, Charles, 
Richardson, and Newitt (ref. 38); Belden and Kassel (ref. 45); and the theses of Hinkle 
(ref. 31); Helander (ref. 35); Mitlin (ref. 36); and Depew (ref. 39). The experimental 
data from these works were thoroughly analyzed by Julian (ref. 46) in order to verify his 
eddy viscosity model. Some data were corrected for acceleration or static effects, and 
other inconsistent data were rejected (e. g. , Belden and Kassel's data were rejected in 
their entirety). Therefore, a description of these works will not be repeated herein; 
however, a summary of the range of the operating variables involved appears in table m. 
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Figure 17. - Fanning friction factor of suspension as function of gas Reynolds 
number as predicted by reference 32. 


Correlations 

The object of this section is to obtain a working correlation that is descriptive of 
the available data and possesses the simplest form for ease of application. The previous 
data sources have indicated that the most significant factors which should be considered 
in a correlation for pressure drop of a gas-solid suspension are loading ratio r], gas 
Reynolds number Reg, and particle diameter D p . The application of the drag coeffi- 
cient model not only appears difficult since V p is not a readily available quantity but 
leads to incorrect results as illustrated by the work of Peskin and Dwyer; and Mehta, 
Smith, and Comings. Therefore, the subsequent analysis will be directed entirely to- 
wards the equivalent friction factor model in order to obtain a working correlation for 
fg/fg, which then can be used to predict suspension pressure drops. 

Equivalent frictio n factor based on eddy viscosity model. - In the eddy viscosity 

model proposed by Julian and Dukler (ref. 32), the value of the friction factor f' 

^ „ s 

depends directly on the value of K, the modified von Karman constant that appears in 
their expression for the eddy viscosity. On the basis of the experimental data of six 
independent investigations of gas -solid suspension pressure drop, Julian and Dukler con- 
cluded that K is a function of loading ratio rj only (see eqs. (61) and (62)) and is inde- 
pendent of gas Reynolds number, particle diameter, and particle density. Based on the 
value of K given by equation (61), the dependence of f' on gas Reynolds number and 
loading ratio was obtained by Julian (ref. 32) and is plotted in figure 17. The plot shows 
that fg increases with solid loading ratio and decreases exponentially with increasing 
gas Reynolds number at constant rj. The slope of the lines of constant tj is approxi- 
mately equal to -0.2, which indicates that the effect of gas Reynolds number on the sus- 
pension flow is approximately the same as it is for pure gas flow (eq. (45)). Thus, for 
the range of conditions covered, fg/f g is essentially independent of Reg. 

A cross plot of figure 17 giving f^ as a function of 1 + fj with Re g as a parameter 
is presented in figure 18. For simplicity, these curves can be fitted approximately by 
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Figure 18. - Fanning friction factor of suspension as func- 
tion of loading ratio plus one (17 + 1 ) as predicted by ref- 
erence 32. Dashed lines are plots of equation (69a). 


the equation 

fg = 0. 046 Re g °' 2 (1 + 77 )°* 3 (69a) 

denoted by the dashed lines in figure 18. 
Thus, a very simple approximation for the 
Julian and Dukler correlation can be ob- 

T . 

tained in terms of f /f : 
s' g 

f’ 

-® = (1 + if) 0 * 3 (69b) 

f 


Equivalent friction factor by use of Reynolds analogy . - The use of the Reynolds 
analogy in conjunction with the correlations of the heat-transfer coefficient represented 
by equations (44) and (43) yields four different friction factor correlations. Applying the 
Reynolds analogy separately to the particles and then superposing the two as in the man- 
ner of Gorbis (eq. (55)) yields the following two correlations: 


— = 1 + 4.0 Re" 0, 32 n 
f S ' 

g 


(70) 


and 


— = 1 + ± 
f 6 

g 


7. 6 ( 677 )°* 45 Re~°‘ 21 - 1 

B 


(71) 


Equation (70) is based on equation (44) for h/h , whereas equation (71) is based on 

s & 

equation (43) for h s /h g . In a similar manner, when the Reynolds analogy is applied di- 
rectly to the suspension as given by equation (57), the correlations obtained from equa- 
tions (44) and (43) are, respectively, 


f s_ ( 1 + 4 - 0 Re g 0 - 3 2 6 ?? ) (1 + r\) 

t g 1 + 


P2) 


and 
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Figure 19. - Comparison of Reynolds analogy correlations. Reynolds number, 25 000. 


fg_ 7.6(^) 0 ' 45 Re”°‘ 21 (1 + rj) 


(73) 


The four correlations obtained (eqs. (70), (71), (72), and (73)) are plotted as a func- 
tion of ?7 with 6 as a parameter in figure 19 for a gas Reynolds number of 25 000. 
Figure 19(a) shows the correlations for the superposition model. With the superposition 
model and the heat- transfer correlation that leads to equation (70), the ratio of the 
equivalent friction factor of the suspension to the friction factor of the gas is independent 
of the specific heat ratio 6. Furthermore, when the heat-transfer correlation associ- 
ated with equation (71) is used, the dependence of the equivalent friction factor ratio on 
6 is not large, especially for loading ratios between 4 and 10. 

Figure 19(b) shows that the suspension model yields a considerable dependence of 
the equivalent friction factor ratio on the specific heat ratio. Use of either heat-transfer 
correlation gives about the same results. It should be pointed out that fg/f g for equa- 
tion (72) with 6 = 2 drops below unity for values of 77 < 2 before rising again to unity 
at 77 = 0. 

A comparison of figures 19(a) and (b) shows that the two methods of applying the 
Reynolds analogy yield exactly the same results when 6 = 1; that is, equation (72) 
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Figure 20. - Comparison of superposition Reynolds analogy with eddy viscosity model. 


reduces to equation (70), and equation (73) reduces to equation (71). 

Of the four possible correlations for f' /f , equation (70) is clearly the simplest, 

® D 

since it does not show a dependence on 6. It does not seem probable that the specific 
heats of the particles and the gas should affect pressure drop; in fact, specific heats 
have not been considered as a factor in any of the pressure drop literature. Therefore, 
until a definite effect of 6 on pressure drop is established, it seems reasonable to con- 
sider only equation (70) to represent the Reynolds analogy. Equation (70) is also in the 
form of the original Gasterstadt equation (eq. (32)), which has been used for correlating 
pressure drop data by a number of investigators (refs. 42, 43, and 44) and which provides 
a ratio of unity at rj = 0. 

When equation (70) is multiplied by f (eq. (45)), 

O 


f’ = 0. 046 Re *°- 2 + 0. 184 Re‘°* 52 v 
® & & 


(74) 


Figures 20(a) and (b) compare f g (as calculated by eq. (74)) to the curves obtained by 
Julian and Dukler (fig. 17). As can be seen from figure 20(a), the friction factor as 
predicted by the Reynolds analogy correlation shows a greater dependence on Re g as r\ 
is increased than that found by Julian and Dukler. 

At very low values for the loading ratios the first term on the right-hand side of 
equation (74) predominates, and f g is proportional to Re”^' 2 . At very high values for 
the loading ratios the second term predominates, and f’ is proportional to Re - ®' ® 2 . 
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Figure 21. - Ratio of Fanning friction factor of suspension to friction factor of gas as 
function of loading ratio for various correlations and data sources. 

Therefore, the Re^ dependence of f^ varies between the -0.2 and -0. 52 power of Re^. 

Julian and Dukler found the dependence to remain at approximately the -0. 2 power, as in 

the case of a pure gas. Nevertheless, the agreement between the eddy viscosity model 

and the Reynolds analogy is quite close at Reg > 25 000. 

It should be noted that most of the data upon which Julian and Dukler' s correlation 

is based utilized gas Reynolds numbers ranging from approximately 30 000 to 100 000. 

The only data where lower Reynolds numbers were employed were the data of Depew and 

that of Hariu and Molstad. The data of Depew (Re = 13 500 to 27 400) yield consistently 

higher results for f' than those predicted by Julian and Dukler. The data of Hariu and 
s 

Molstad (Re= 3700 to 11 000) corrected for acceleration effects by Julian, however, 
appear to be in agreement with the correlation. 

Comp arison of models with experimental data . - Figure 21 is a plot of f^/f against 
7j as predicted by the Reynolds analogy (eq. (70)) with gas Reynolds number as a param- 
eter. Julian and Dukler's eddy viscosity model appears as a single (dashed) curve in- 
dependent of gas Reynolds number. The data of Mitlin; Hinkle; Clark, Charles, 
Richardson and Newitt; Depew; Hariu and Molstad; and Helander are adequately represen- 
ted by this curve. The data of Mehta, Smith, and Comings (36 -m particles), Vogt and White's 
clover data, and points from Farbar's curve (see fig. 16) are also presented in the figure. 
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Figure 22. - Fanning friction factor of suspension as function of particle diameter. Load- 
ing ratio, 3 < 7) < 4. 


The agreement of these data with both the Reynolds analogy and the eddy viscosity 
model curve is very close, although a Reynolds number effect similar to that predicted 
by the Reynolds analogy seems to be discernible; that is, the data points at low gas 
Reynolds numbers are somewhat higher than those at high gas Reynolds numbers. These 
results are also in agreement with the correlation presented by Dogin and Lebedev 
(see eq. (67b)), where a similar type of gas velocity dependence is predicted. 

At values of rj between 4 and 10, figure 21 shows that the friction factor f* /f 

s g 

curves may be fitted by straight lines of slope of the order of 0. 45 as was found for the 
heat-transfer coefficient correlation (eq. (43)). The data of Peskin; Mehta, Smith, and 
Comings; and Farbar also show approximately this same r\ dependence on the friction 
factor (see figs. 14 to 16), which is a good indication of the validity of the analogy. 

Peskin’ s data; and Mehta, Smith, and Comings’ 97-micron particle data have not 
been plotted in figure 21 since the absolute magnitude of their data is much higher than 
predicted by the eddy viscosity model or the Reynolds analogy and is in disagreement 
with most of the other data previously discussed. 

Effect of particle diameter . - Since the heat-transfer correlations utilized (eqs. (43) 
and (44)) showed essentially no particle diameter dependence, the friction factor cor- 
relations naturally show no particle diameter effect. This fact is in agreement with the 
findings of Julian and Dukler, who could discover no apparent particle diameter effect 

on f for the data that they investigated, which included a particle diameter range from 

® 1 

25 to 6000 microns. Figure 22 is a plot of f against particle diameter obtained from 
the experimental data for loading ratios between 3 and 4. Although the data do not fall 
on a single horizontal line (which would indicate no particle diameter effect), no other 
obvious particle diameter effect is evident. Therefore, the conclusion that particle 
diameter has no effect as predicted by both the eddy viscosity model and the Reynolds 
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analogy must be accepted until further data become available. 

Final correlations . - As a result of the preceding analysis, it would appear that for 
the range of variables, covered in the literature either the following correlation obtained 
using the Reynolds analogy 


f’ 


-= 1 + 4.0Re'°' 32 ?7 

& 


g 


(70) 


or the eddy viscosity model as proposed by Julian and Dukler, which can be expressed 
approximately as 



f 


g 


(1 + rj) 


0.3 


(69b) 


yields values of f’ /f , which are in agreement with most of the experimental data. Both 
s g 

correlations are therefore recommended for predicting pressure drops associated with 
the fully developed turbulent flow of gas-solid suspensions in smooth tubes. 


APPLICATION OF RESULTS 

The correlations obtained as a result of this analysis should be applied in flow system 
calculations with heat transfer only after first examining the various limitations implicit 
in their development. Aside from the restrictions with respect to gas Reynolds number, 
particle diameter, and loading ratio previously indicated there are other less obvious 
restraints and implications. These factors arise from the fact that the experimental data 
from which the correlations were developed were obtained for certain combinations of 
gases and particles at certain test conditions; that is, fixed pressures, temperatures, 
heat fluxes, velocities, and pipe diameters. 

For the heat-transfer coefficient, the test data in the survey covered pressures 
from 0 to 130 pounds per square inch gage, temperatures from 75° to 1100° F, velocities 
from about 1 to 200 feet per second, and pipe diameters from 0. 3 to 4 inches. In practi- 
cally all cases, heat was added to the suspension flow. For the friction factor, all tests 
were run at essentially adiabatic conditions at low temperatures (~70° F). The range of 
pressures, velocities, and pipe inside diameters for this factor were, respectively, 

0 to 30 pounds per square inch gage, 10 to 150 feet per second, and 0. 25 to 5 inches. 

For both the heat-transfer and friction factor tests, the flow could be considered as fully 
developed turbulent flow. 
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In applying the correlations obtained herein to calculations for any other system the 
following procedure is recommended. The following equation may be used for the sus- 
pension heat-transfer coefficient 

h ^(r) \ ' (75) 

\ g/corr 


where (h /h ) is obtained from the appropriate correlation presented herein, 
s g corr 

Thus, all that need be done is to compute ( h s / h g) corr and both evaluated at the 
temperature, pressure, velocity, and pipe diameter associated with the gas in the sus- 
pension flow. It should be pointed out, however, that in various application systems all 
these variables may not be independent. For example, in a Brayton power cycle, the 
local pressure, temperature, and volumetric flow rate are interrelated through the cycle 
input parameters. Furthermore, these local state values may in turn also be functions 
of the suspension loading ratio. 

In a manner similar to that used for the suspension heat-transfer coefficient, the 
suspension effective friction factor can be found from 



(76) 


where both f and (f* /f ) are to be evaluated in the same manner as outlined for 
g s' g corr 

(h /h ) and h . It is apparent, therefore, that in using the preceding procedure, 
s g corr g 

correlations of the form of equations (44) and (70) are preferred as they conveniently 
reduce to unity when 77 = 0. 

Strictly speaking, in equations (75) and (76), both h g and f refer to the values of 

h and f’ for 77 = 0; however, for practical purposes, h and f can be evaluated 
s s & & 

from appropriate, available pure gas correlations. 

The question also exists of the applicability of the selected correlations for 

(h/h )„ and(f’/f L to gases, particles, state conditions, and heat-transfer 

' s' g corr s' g corr 

methods different from those covered in the referenced experimental data. However, 
in view of the approximate nature of the correlation and the large scatter in the refer- 
enced data, such questions may not be serious for preliminary calculations. 

In order to increase the correlation precision, further work over a wide range of 
gases, particles, heat-transfer methods, heat fluxes, temperatures, pressures, veloc- 
ities, and pipe diameters is desirable. This experimental work should include pure gas 
(77 = 0) data with the same apparatus because these data may not be in exact agreement 
with the established pure gas correlations. 
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SUMMARY OF RESULTS 




The published literature data concerning the heat-transfer and pressure drop 
associated with the flow of dilute gas-solid suspensions were analyzed. Although the 
available data were insufficient and somewhat inconsistent for precise correlation, ap- 
proximate relations for preliminary use in predicting suspension heat-transfer coeffi- 
cients and friction factors were developed. Two possible correlations were developed 
for the heat-transfer coefficient ratio. From these correlations and the available data, 
it was found that 

1. The ratio of the heat-transfer coefficient of a suspension to the heat- transfer 
coefficient of an equal volumetric rate of flow of gas at the same pressure and temper- 
ature increases with loading ratio (lb solid/lb gas). 

2. The heat-transfer coefficient ratio increases as the ratio of particle specific heat 
to pure gas specific heat is increased. 

3. The heat-transfer coefficient ratio decreases as the gas Reynolds number is in- 
creased. 

4. The particle diameter has little effect on the ratio of heat-transfer coefficients 
for the range of particle sizes between 30 and 150 microns. 

Application of the Reynolds analogy to these heat-transfer correlations yielded a 
relation for the ratio of equivalent friction factor of a suspension to the equivalent fric- 
tion factor of a gas. This relation indicates the following: 

1. The friction factor ratio increases as the loading ratio is increased. 

2. The friction factor ratio decreases as the gas Reynolds number is increased. 

3. The particle diameter and specific heat appeared to have little or no effect on 
the friction factor. 

The results of experimental data available in the literature as well as an analytical 
study of the equivalent friction factor using an eddy viscosity model for the suspension 
compare quite favorably with the developed friction factor correlation especially at low 
loading ratios (77 < 5). However, it was evident from the correlations of both the heat- 
transfer coefficient and the equivalent friction factor data that additional experimental 
investigations may be required to obtain more precise variations for specific particles 
and variable ranges of interest. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 12, 1966. 
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APPENDIX - SYMBOLS 


a 

» 

a 

B 

b 

C 



F 
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constant (eq. (35)) 

experimentally determined con- 
stant (eq. (63)) 

constant (eq. (68)) depending 
upon solid-gas system 

constant (eq. (36)) 

constant (eq. (67b)) depending 
upon gas -solid system 

standard drag coefficient for 
spherical particles 

specific heat of gas, 

Btu/(lb)(°F) 

specific heat of particles, 
Btu/(lb)(°F) 

specific heat of suspension, 
Btu/(lb)(°F) 

tube inside diameter, ft 

diameter of particle, p. or ft 

complicated function of gas 
Reynolds number and other 
parameters in Gasterstadt 
equation (eq. (32)) 

drag force, lb 

Froude number for gas 

Fanning friction factor for pure 
gas 

mixture Fanning friction factor 
(eq. (63)) 

particle Fanning friction factor 
based on velocity of particle 
and density of gas 


fp particle Fanning friction factor 

based on velocity of particle 
and particle density in mix- 
ture 

f Fanning friction factor for sus- 

s 

pension based on velocity and 
density of suspension 

fg Fanning friction factor based 

on velocity of suspension and 
density of gas 

2 

g acceleration of gravity, ft/ sec 

g c gravitational constant, 

(lb mass)(ft)/(lb force)(sec 2 ) 

h„ convective heat-transfer coef- 

g 

ficient for pure gas, 
Btu/(hr)(ft 2 )(°F) 

hp convective heat-transfer coef- 

ficient, Btu/(hr)(ft 2 )(°F) 

(eq. (27)) 

h convective heat-transfer coef- 

s 

ficient for suspension, 

Btu/ (hr ) (f t 2 ) (° F) 

K von Karman constant for a pure 

gas and a function of 77 for a 
suspension 

k experimentally determined con- 

stant (eq. (60)) 

kg thermal conductivity of pure 

gas, Btu/(hr)(ft)(°F) 

k thermal conductivity of particle, 

Btu/ (hr) (ft) (°F) 

k g thermal conductivity of sus- 

pension, Btu/(hr)(ft)(°F) 


1 11 mi mu 1 mm 
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k i’ k 2>l 

constants when all parameters 

T 

m 

mean temperature of wall and 

k 3 ,k 4j 

except r] are fixed 

suspension, °R 

k l, k 2’l 

constants when all parameters 

T 

n 

temperature of particle, °F 

k 4 J 

' except 6 and r\ are fixed 

r 

T s 

bulk temperature of suspension, 

L 

length of tube, length of 

°R 

£ 

calming plus test sections, 
ft 

V g 

V n 

gas velocity, ft/sec 
velocity of particle, ft/sec 

exponent in equation (68) and 

P 


function of dimensionless 
1/2 

group, Re p (C d ) 

V PS 

terminal settling velocity of 
particle, ft/sec 

m 

experimentally determined 

v s 

velocity of suspension, ft/sec 


constant (eq. (60)) 

w„ 

g 

mass flow rate of pure gas. 

Nu_ 

g 

Nusselt number of pure gas 

lb/sec 

Nu 

s 

n 

Nusselt number of suspension 
number of particles per cubic 

W 

P 

mass flow rate of solids, 
lb/sec 


foot of mixture, 
particles/ft^ 

y 

radial distance from tube wall, 
ft 

AP 

g 

ap p 

pressure drop caused by gas 
flow, lb force/ft^ 

additional pressure drop 

4 

dimensionless distance from 
tube wall evaluated at center 
of tube 


caused by adding particles 

2 

to gas stream, lb force/ft 

a 

g 

sonic velocity of pure gas, 
ft/sec 


pressure drop caused by 
suspension flow, 

a 

s 

sonic velocity of suspension, 
ft/ sec 

[ 1 + ( VV]/[ l -VV] 


lb force/ft^ 

0 

Pr g 

Prandtl number of pure gas 

r g 

isentropic specific heat ratio 

Pr s 

Prandtl number of suspension 

for pure gas 


Reynolds number of pure gas 

^s 

isentropic specific heat ratio 

Re p 

particle Reynolds number 


for suspension 

Re s 

Reynolds number of 

6 

ratio of specific heat of particle 

suspension 


to specific heat of gas 

T 

g 

temperature of gas, °F 

e 

fractional solid volume, 
ft^ solid/ft^ suspension 
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n 

e' eddy viscosity, ft /hr 

T) loading ratio, 

lb mass solid/lb mass gas 

H viscosity of pure gas, 
lb/ (ft) (hr) 

jU g viscosity of suspension, 
lb/ (ft) (hr) 

Pg density of pure gas, 
lb/ft 3 gas 


O 

p particle density, lb/ft solid 

Jr 

density of solid, 
lb solid/ft 3 suspension 

p_ bulk density of suspension, 

S Q 

lb/ft suspension 

(p constant, (y + m - 60)/22 

xjy factor characteristic of loading 

ratio (eq. (7a)) 
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